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Abstract 

In the context of the currently ongoing efforts to improve the ac- 
curacy and reliability of the measurement of the Lense-Thirring effect 
in the gravitational field of the Earth it has recently been proposed 
to use the data from the existing spacecraft endowed with some active 
mechanisms of compensation of the non-gravitational accelerations like 
GRACE. In this paper we critically discuss this interesting possibility. 
Unfortunately, it turns out to be unpracticable because of the impact 
of the uncancelled even zonal harmonic coefficients of the multipolar 
expansion of the terrestrial gravitational potential and of some time- 
dependent tidal perturbations which would resemble as superimposed 
linear trends over the necessarily limited observational time span of 
the analysis. 

1 Introduction 

In this paper we analyze a recent proposal of using the data from the 
existing spacecraft with some active mechanisms of compensation of the 
non-gravitational perturbations to improve the ongoing tests concerning the 
general relativistic Lense-Thirring effect in the gravitational field of the 
Earth. 



2 The proposal of using the existing satellites en- 
dowed with accelerometers 

At the end of Section 5, p. 649, right column of it is written: "[...] one 
wonders why the author of Iorio 1 (2005) has not proposed the use of the 
1 It is the reference [3] of this paper. 
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GRACE satellites which not only have by far better shape and orbital sta- 
bility compared to JASON but they also carry ultra precise accelerometers 
that measure all non-conservative accelerations, so that they are in practice 
"free-falling" particles in vacuum, at least to the extent that is covered by 
the accuracy of these instruments." 

The answer is that the existing spacecraft with some active mechanism of 
compensation of the non-gravitational forces (CHAMP, GRACE, GP-B) fly 
at too low altitudes in polar orbits, so that the use of their nodes, combined 
or not with those of the LAGEOS satellites, would greatly enhance the sys- 
tematic errors both due to the uncancelled even zonal harmonic coefficients 
Ji of the multipolar expansion of the terrestrial gravitational potential and 
to certain time-dependent tidal perturbations. 

2.1 The use of the node of GP-B 

In fact, the possible use of the node of the drag-free spacecraft GP-B has 
been considered in 2 [3]. The major problems come, in this case, from the 
fact that a nearly perfect polar orbital configuration would make the node to 
precess at a very low rate. This is a problem because one of the major time- 
dependent tidal perturbation, i.e. the solar K\ tide, which is not cancelled 
out by the linear combination approach because it is a tesseral (m = 1) 
perturbation [5], has just the same period of the satellite's node. So, over an 
observational time span of about one year-which is the lifetime of GP-B-it 
would resemble an aliasing superimposed secular trend 7 . For example, in 
[2] a combination with the node of LAGEOS and GP-B has been considered: 
the coefficient of GP-B would amount to -398, thus fatally enhancing such 
tidal bias. 

2.2 The use of the nodes of the other existing geodetic satel- 
lites and of CHAMP and GRACE 

Analogous conclusions can also be drawn for a possible use of CHAMP 
and GRACE. Indeed, let us consider, e.g., a combination with the nodes 
of LAGEOS, LAGEOS II, Ajisai, Starlette, Stella, CHAMP and one of the 
two spacecraft of the GRACE mission in order to cancel out the first six 
even zonal harmonics Ji,£ = 2,4,6,8, 10, 12. Apart from the coefficient of 
LAGEOS, which is 1, the coefficients which weigh the nodes of the other 

2 Another analysis can be found in 
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The slope of the gravitomagnetic trend is 3714.58 milliarcseconds per year 
(mas yr" 1 in the following). The systematic error due to the uncancelled 
even zonal harmonics J14, Jiq, Jis, ... amounts to 4.4%, at 1-sigma level, ac- 
cording to a calculation up to degree t = 42 with the variance matrix of the 
Earth gravity model EIGEN-CG01C [S] which combines data from CHAMP, 
GRACE and ground-based measurements. Such error level is not compet- 
itive with those which can be reached by the combination with LAGEOS 
and LAGEOS II proposed in 3 |101 1111 112j and used by Ciufolini, Pavlis and 
Peron in their recent test and the combination involving also Ajisai 
and Jason-1 |141 116j . Indeed, EIGEN-CG01C yields for them a systematic 
error of gravitational origin of about 6% and less than 2%, respectively, at 
1-sigma level. Moreover, while the error due to the geopotential is sensitive 
to the even zonals, at most, up to degree I = 20, for such combinations, the 
multi-combination with CHAMP and GRACE is sensitive to a much larger 
number of even zonal harmonics due to the inclusion of the lower altitude 
satellites Starlette and, especially, Stella, CHAMP and GRACE. This makes 
rather difficult and unreliable the evaluation of the systematic bias induced 
by the static part of the geopotential because the calculation of the coeffi- 
cients of the classical secular precessions |17j yield unstable results after 
degree I ~ 40. Thus, the 4.4% estimate of the bias due to the even zonals is 
probably optimistic. Another serious drawback of such multi-combination 
is represented by the relatively large magnitude of the coefficients 
which weigh the nodes of CHAMP and GRACE. Indeed, they enhance the 
impact of all the uncancelled time-dependent perturbations among which 
the solar K\ tide is one of the most powerful. The periods of the related 
orbital perturbations amount to -2.63 years for CHAMP and -7.23 years for 
GRACE. Such effects would represent serious aliasing bias over the neces- 

3 The idea of using only the nodes of LAGEOS and LAGEOS II in view of the expected 
improvements of our knowledge of the terrestrial gravitational field with GRACE was put 
forth for the first time by Ries et al. in 9; . 
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sarily short observational time span due to the limited lifetimes of CHAMP 
and GRACE with respect to the geodetic satellites. 

2.3 The use of the nodes of the other existing geodetic satel- 
lites 

In regard to the possibility of only using the other existing geodetic spherical 
satellites, mainly Ajisai, Stella and Starlette due to their long data records 
available, this problem has already been tackled in a number of papers, 
like |181 I19j and [2U|. The fact that the inclusion of the other satellites 
in the linear combination scheme is still not competitive, although the im- 
provements due to the first models from CHAMP and GRACE, was shown 
in [Tl]: the systematic error of gravitational origin amounts to 31% with 
EIGEN-CG01C. 

2.4 The impact of the new Earth gravity models from CHAMP 
and GRACE on the use of a polar orbital geometry 

The impact of the most recent Earth gravity models from CHAMP and 
GRACE on the use of the node of a single satellite in polar orbit has been 
discussed in detail in ^Sj. By using EIGEN-CG01C, it turns out that for 
a semimajor axis of, e.g., 8000 km, quite larger than GP-B, CHAMP and 
GRACE, the systematic error due to the full spectrum of the even zonal 
harmonics would amount to 25% for an inclination of 88 deg. In this case 
the period of the K\ tide would amount to ~ 10 3 days. Instead, for an 
inclination of 89.9 deg the bias due to the even zonals would be 2% but the 
period of the tidal perturbation would raise to ~ 10 days. Also with the 
new terrestrial gravity field solutions the linear combination approach would 
fail. 

3 The proposed use of Jason- 1 

Iorio proposed to investigate the possibility of analyzing a suitable linear 
combination involving the nodes of LAGEOS, LAGEOS II, Ajisai and Jason- 
1 in 4 HU. 

The advantages of the combination involving also Jason-1 are 

• Cancellation of the first three even zonal harmonics J2, J4, Jq of the 
geopotential along with their secular variations J2, J4, J§- Moreover, 

4 See also on the same topic. 
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the systematic error due to the remaining higher degree even zonal 
harmonics Jg, Jio, ••• is almost model-independent: indeed, it is < 2% 
(1-sigma level) according to the 2nd generation GRACE-only Earth 
gravity models EIGEN-GRACE02S |H] (2%) and GGM02S (2.7%) 
and to the model EIGEN-CG01C 1.6%). It should be mentioned 
that it is expected that GRACE will yield a larger improvement in 
the knowledge of the higher degree even zonal harmonics, to which 
the Jason's combination is sensitive, instead of the lower degree even 
zonals, which mainly affect the node-node LAGEOS-LAGEOS II com- 
bination. This means that it may happen that, in the near future, the 
bias due to the geopotential will reduce down to ~ 1% for the Jason's 
combination in a satisfactorily model-independent way, while it may 
remain more or less unchanged for the two-nodes LAGEOS-LAGEOS 
II combination. The latest results obtained with the combined model 
EIGEN-CG03C [22] seems to confirm this trend [23]. 

Finally, only the first ten even zonal harmonics (i.e. up to J20) should 
be accounted for in the sense that the systematic error due to the 
geopotential does not change with the inclusion of more zonals in the 
calculation. 

• Small coefficient, ~ 1CP 2 , of the node of Jason-1. This is particularly 
important for reducing the impact of the non-gravitational acceleration 
suffered by Jason-1. 

• No secular or long-period perturbations of gravitational and non-gravitational 
origin. 

The possible weak points of this proposals are mainly the following 

• The huge impact of the non-gravitational forces-mainly atmospheric 
drag, direct solar radiation pressure and Earth's albedo-on Jason-1 
which has not a spherical shape, being endowed with solar panels. 
Moreover, they should be modelled in a truly dynamical way in order 
to avoid to absorb the Lense-Thirring effect as it would happen in the 
empirical reduced-dynamic approach adopted so far. 

• The difficulty of getting a smooth long time series of its orbit also due 
to the periodical orbital maneuvers. On the other hand, it should be 
noted that, due to the Jason's main goal which is ocean altimetry, 
only the radial and along-track in-plane components of its orbit have 
received major attention up to now. Also the orbital maneuvers are 
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mainly, although not entirely, in the orbital plane. Instead, the node 
is related to the cross-track, out-of-plane component of the orbit. 

However, a detailed, although preliminary, evaluation of the impact of the 
Jason's non-conservative forces on the entire combination has been per- 
formed in |14| . No secular effects occur. On the contrary, a relatively high- 
frequency (the 120-days period of the (3 cycle related to the orientation 
of the orbital plane with respect to the Sun) of non-gravitational origin has 
been found. Its impact on the suggested measurement of the Lense-Thirring 
effect has been evaluated to be < 4% over a 2-years time span (without re- 
moving such periodic signal). 

4 On some possible misunderstandings concern- 
ing various aspects of the measurement of the 
Lense-Thirring effect with the LAGEOS satel- 
lites. 

Among the six Keplerian orbital elements in terms of which it is possible 
to parameterize the orbital motion of a test particle in the gravitational 
field of a central body of mass M [231, the longitude of the ascending node 
Q, the argument of the pericentre u and the mean anomaly A4 undergo 
secular precessions due to the even zonal harmonics. The general relativistic 
gravitomagnetic force only affects Q and uj with the secular precessions of the 
Lense-Thirring effect In principle, Q,uj and Ai could all be used in order 
to design suitable linear combinations in order to cancel out the low-degree 
even zonal harmonics whose classical precessions are much larger than the 
Lense-Thirring signals of interest. However, the mean anomaly cannot be 
used at all because it is sensitive to huge non-gravitational perturbations 
which especially affect the Keplerian mean motion n = w GM/a 3 via the 
indirect effects on the semi-major axis a. They have quadratic signature in 
time. 

The title of Section 5, p. 648 of is: "On the use of the mean anomaly 
and on the use of Jason to measure the Lense-Thirring effect proposed in 
Iorio (2004) ". Moreover, at the beginning of Section 5, p. 648, right column, 
first paragraph of 1 it is written that "[...] one of the most profound 
mistakes and misunderstandings of Iorio 5 (2005) is the proposed use of the 
mean anomaly of a satellite to measure the Lense-Thirring effect (in some 

5 It is the reference of this paper. 
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previous paper by the same author the use of the mean anomaly was also 
explicitly proposed)." 

In regard to the mean anomaly, the work quoted by the authors of 
in the title of Section 5 of their paper is an old, unpublished version of the 
preprint gr-qc/0412057vl. When the paper was submitted (23 March 
2005) and resubmitted (4 April 2005) the version v3 of gr-qc/0412057 , with 
substantial revisions, was already available on the Internet since 22 Febru- 
ary 2005. In all the versions of such a preprint the mean anomaly is not 
mentioned at all. Moreover, no papers by Iorio in which such an alleged 
explicit proposal would appear are quoted in the text of Section 5 of [Q. In 
fact, in the whole literature does not exist any published paper by 6 Iorio 
in which he explicitly put forth the possibility of using the mean anomaly 
of the LAGEOS satellites for measuring the Lense-Thirring effect. On the 
contrary, the use of the node and the perigee of the Earth' satellites is always 
explicitly mentioned in many papers where the linear approach combination 
is exposed and generalized to other situations (see, e.g. pp. 16-17 of [T8] . 
pp. 5475-5476 of |H], p.3 of US], pp. 289-290 of [T7j, p. 5 and p. 8 of Q3], 
pp. 4-5 of [__0])- Moreover, in Section 2.1.2, p. 607, end of the left column 
of [Sj , after a general description of the linear combination approach for the 
measurement of the Lense-Thirring effect, it is unambiguously written "In 
general, the orbital elements employed are the nodes and the perigees and 
the even zonal harmonics cancelled are the first N-l low-degree ones." In 
Section 5, p. 649, left column, second paragraph of 4_ it is written: "A more 
detailed work discussing [...] the use of the mean anomaly and other highly 
unfeasible proposals by Iorio (2005) will be the subjects of following paper" . 
It is, thus, likely that the authors of pQ, after a careful inspection of the 
relevant literature, will experience some difficulties in finding material for 
such an announced paper, at least as far as the use of the mean anomaly for 
measuring the Lense-Thirring effect is concerned. 

Instead, it turns out that one of the authors of [Q used the semimajor 
axis, the eccentricity and the mean anomaly of LAGEOS II in some (not 
explicitly released) linear combination with the nodes of LAGEOS and LA- 
GEOS II and the perigee of LAGEOS II in some previous tests with the 
EGM96 Earth gravity model (Figure 4 and pp. 47-48 of [27]; Figure 4 and 
paragraph (b), pp. 2376-2377 of |5B]). 

In Section 4, p. 647 of it is claimed that one of the authors would have 
put forth the idea of suitably combining the Keplerian orbital elements of 

6 The updated list of the accepted and published papers by Iorio is available on the 
Internet at http:/ /digilander. libero.it/lorri/list_of_publications. htm 
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the LAGEOS-type satellites in order to disentangle the Lense-Thirring effect 
from some of the much larger biasing secular precessions induced by the even 
zonal harmonics of the geopotential since the beginning of his studies on the 
measurability of the gravitomagnetic effect, when the LAGEOS II, launched 
in 1992, did not exist yet. This claim is supported by quoting from [2l"l 13*0] : 
"[...] A solution would be to orbit several high-altitude, laser-ranged satel- 
lites, similar to LAGEOS, to measure J2, J4, Jq, etc., and one satellite to 
measure ^Lense-Thirring-" It seems more likely that this statement could, at 
most, mean that one should first use many satellites to accurately measure 
the various Earth's even zonal harmonics with a sufficiently high accuracy 
and, then, analyze the node only of one satellite for safely measuring its 
gravitomagnetic Lense-Thirring precession. However, there is no trace at all 
of the linear combination approach which will be introduced only in 1996 
|31| in the particular case of the nodes of LAGEOS and LAGEOS II and the 
perigee of LAGEOS II. Indeed, the previously quoted statement comes after 
a discussion of the impact of J2 and of the higher degree even zonal har- 
monics on the possible use of the node only of LAGEOS; moreover, the rest 
of the papers |29| I30| deals with the supplementary configuration LAGEOS- 
LAGEOS X (later LAGEOS-LAGEOS III /LARES /WEBER-SAT) . 

The remarks about the triviality of the simple algebraic linear system 
of two equations with two unknowns with which the LAGEOS-LAGEOS II 
combination is obtained seem to be a little inappropriate because they could 
also be extended to the original node-node-perigee combination and the 
related linear algebraic system of three equations with three unknowns. 

In regard to the combination of the nodes of the LAGEOS satellites, 
the papers on it are available on the Internet since April 2003. Moreover, 
their author, who personally knows the authors of ""Q having collaborated 
with them for some years, sent them various e-mails 7 between April and 
September 2003 with his results attached. In one case (28 March 2003 
and 30 March 2003), Ciufolini asked Iorio to prepare a short table with his 
results in .doc format in view of a video-conference with NASA'officials to 
be attended in the following days by Ciufolini. A few days before, 26 March 
2003, Iorio e-mailed the .pdf of 10 to Ciufolini. Later, on 7 September 
2003, Iorio discussed the impact of the GGM01C Earth gravity model on 
the LAGEOS-LAGEOS II combination with Pavlis. 

Thus, we feel that the style and the content of the comment in the middle 
of left column, p. 648 of ""Q: "In conclusion, all the claims of Iorio (2005) 
are simply lacking of any rational basis: above is shown how much work was 

7 They are all available on request. 
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already published on this topic before Iorio (2005) even began to produce 
any of his paper on this topic and to rediscover some earlier results [...]. 
To avoid the misunderstandings of Iorio (2005), it would have just been 
a matter of very careful reading the previously existing literature on this 
subject!" sound rather inappropriate. 

5 Conclusions 

In this paper we demonstrated that the interesting proposal of using the or- 
bital data from the existing terrestrial spacecraft endowed with some active 
mechanism of compensation of the non-gravitational perturbations like GP- 
B, CHAMP and GRACE is not suitable for measuring the Lense-Thirring 
effect both because of their too low altitude and of their nearly polar orbital 
geometry. Indeed, they would introduce in the resulting linear combina- 
tions much more even zonal harmonics of high degree (more than t = 40) 
which would enhance the induced systematic error and make unreliable its 
calculation. Moreover, the polar geometry of GP-B, CHAMP and GRACE 
would have, as a consequence, that the coefficients with which they would 
enter the combinations would be larger than 1, thus enhancing various un- 
cancelled time-dependent long-period perturbations like that due to the K\ 
tide. The new Earth gravity models from CHAMP and GRACE do not alter 
this situation. 
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